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The  photon  correlations  of  resonance  fluorescence  emitted  by  a  degenerate 
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optical  laser  drives  a  low-lying  atomic  transition,  and  the  emitted  fluorescence 
is  detected  by  a  photon  counter.  An  expression  for  the  spontaneous -decay 
operator  is  derived,  and  it  appears  that  the  various  substates  of  the  upper  level 
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detailed  elaboration  of  the  j  ■  0  ■*  J  »  1  transition,  for  which  the  two-photon 
correlation  is  obtained  explicitly  as  a  function  of  the  lifetime  A(h),  depending 
on  the  atom-surface  distance  h,  the  laser  linewidth,  the  detuning  from  resonance, 
and  the  dipole-coupling  strength  with  the  external  field. 
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I.  INTRODUCTION 

Irradiation  of  an  atom  by  an  optical  laser  with  a  frequency  in  close 
resonance  with  an  electronic  transition  frequency  (level  separation)  will  cause 
the  atom  to  absorb  photons  from  the  incident  field.  Subsequent  spontaneous  decay 
is  accompanied  by  the  emission  of  fluorescence  radiation,  which  is  amenable  to 
observation  in  an  experiment.  Detailed  properties  of  the  scattered  radiation, 
like  its  intensity,  polarization,  spectral  profile  and  temporal  photon 
correlations  reflect  the  dynamical  evolution  of  the  atom  in  its  environment.  A 
generic  method  to  obtain  information  about  interatomic  potentials,  for  instance, 
is  by  measuring  the  line  shape  of  the  fluorescence,  which  is  emitted  during  a 
collision.  In  a  similar  way  it  should  be  feasible  to  study  the  behavior  of  an 
atom  in  the  vicinity  of  a  surface,  or  the  alterations  in  its  optical  properties 
due  to  the  presence  of  an  optically-active  substrate. 

Resonance  fluorescence  of  a  free  atom  (gas  phase  or  atomic  beam)  in  a  strong 
external  field  has  been  studied  extensively  during  the  last  two  decades. 

1-3 

Theory  predicted  a  three-line  spectrum  for  a  two-state  atom,  which  was 

4-6 

confirmed  by  experiments.  The  non-classical  nature  of  fluorescence  radiation 

7-9 

was  established  by  the  observation  of  antibunching  and  sub-poissonian 

statistics1^  in  the  temporal  distribution  of  the  photon  emissions.  Also  combined 

time-  and  frequency-resolved  photon  correlations  were  found  to  be  in  excellent 

agreement  with  theory.11  ^  As  stated  in  the  first  paragraph,  the 

characteristics  of  the  fluorescence  are  affected  by  the  atom's  surroundings. 

Most  notably  is  the  change  in  lifetime  of  an  excited  atomic  state  due  to 

17  18 

confining  boundaries  for  the  radiation,  like  in  a  cavity,  near  a  surface,  or 
19  20-24 

in  a  waveguide.  In  a  series  of  experiments  with  molecular  dye  films  on  a 

25 

dielectric  substrate,  the  theoretically-predicted  dependence  of  the  decay 
constants  (inverse  lifetimes)  on  the  distance  of  the  molecule  to  the  surface  was 


found  indeed.  In  this  paper  we  consider  an  atom  near  a  metal  surface  (with 
perfect  conductivity),  which  is  illuminated  by  a  strong  non-monochromatic  laser 
beam,  and  we  study  the  correlations  between  the  emitted  fluorescent  photons. 
Recently,  the  spectral  distribution  of  the  fluorescence  was  calculated*  '  for 
the  same  configuration,  in  the  case  where  the  atom  could  be  accounted  for  by  a 
two-state  model,  and  where  a  semiclassical  description  of  the  emission  process 
was  adopted. 

II.  DRIVEN  ATOM 

An  atom  is  situated  at  r  ■  he  ,  h  >  0,  above  the  xy-plane,  which  separates 
the  vacuum  z  >  0  from  a  perfectly-conducting  substrate  z  <  0.  This  system  is 
irradiated  by  a  laser,  with  an  electric-field  component 

i(k-r-w.  t-$(t)) 

E.(r.t)  -  E  Re  c  e  ‘  "  L  ,  (2.1) 

-l  -  o 

> 

where  Eq  is  the  amplitude,  e  the  normalized  polarization  vector  (ft*  »  1),  k  the 
wave  vector,  perpendicular  to  c  (k*c  ■  0),  and  $(t)  *  stochastically  fluctuating 
phase,  which  gives  rise  to  the  laser  linewldth.  At  the  surface  this  incident 
field  is  reflected,  and  the  sum  field,  evaluated  at  the  position  of  the  atom 
becomes 

-i(«lLt+d(t)) 

EjU)  -  2RoRe  e  IcAcos(hks>  ♦  iC|Sin(hks)J  ,  (2.2) 

where  the  polarization  vector  is  divided  into  perpendicular  and  parallel 
components  with  respect  to  the  xy-plane,  and  k  -  k*e  is  the  z- component  of  the 
wave  vector.  For  solids  with  a  finite  penetrability  (dielectrics),  the  factor  in 
square  brackets  is  different,  but  the  time  dependence  of  this  external  field 


remains  the  same. 
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The  laser  frequency  is  assumed  to  be  tuned  in  close  resonance  with  the 

level  separation  >  0  between  two  (possibly  degenerate)  levels  with 

energy  ftu  (excited)  and  flw  (ground).  If  we  suppress  irrelevant  quantum 
e  g 

numbers,  we  can  write  the  atomic  wave  functions  as  |j  m  >  and  |j  m  >,  with  j  and 

e  e  g  g 

■  the  angular-momentum  quantum  numbers.  In  terms  of  the  projectors 


P  ■  \  |j  m  ><j  m  I 
e  £  |Je  e  Je  e' 

m 

e 

P  -  )  |j  m  ><j  m  I 

g  L  '  g  g  g  g1 


m 


(2.3) 


(2.4) 


onto  the  multiplets,  the  atomic  Hamiltonian  can  be  represented  by 


H  «  P  +  P 
a  e  e  g  g 


(2.5) 


Coupling  of  the  atomic  dipole  p  to  the  external  field  is  governed  by  the 
interaction  Hamiltonian 


-i(ui.  t+$(t)) 

-p-Ej^t)  ■  -  Jhn(h)e  +  H.c. 


(2.6) 


where  we  made  the  rotating-wave  approximation  and  introduced  a  'Rabi  operator' 


0(h)  ■  2(Ko/m)P#ly*cLcoa(hk£)  ♦  ip-C|Sin(hk8))Pg 


(2.7) 


which  generalizes  the  Rabi  frequency  (a  real  number)  for  non-degenerate  atoms  in 
free  space.  Notice  that  8(h)  is  not  Hermitian  and  that  it  depends  explicitly 
on  the  normal  distance  h  of  the  atom  to  the  surface.  A  convenient  concept  will 
turn  out  to  be  the  effective  (non-Hermitian)  surface-dipole  moment  p  (x),  defined 


P#(x)  -  2(cos(x)px  +  i  sin(x)pg) 


(2.8) 
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where  x  is  a  dimensionless  normal-distance  parameter.  Then  the  Rabi  operator 
becomes 

0(h)  -  (E  ,  /*)  Pa(u  (hk  )-c)Po  .  (2.9) 

o  e  -s  z  -  g 

The  equation  of  motion  for  the  density  operator  p(t)  of  the  atom  in  the 
laser  field  and  near  the  surface  is 

i,ldt  *  lHa  ‘  '  il4r(h)p  *  (2.10) 

where  the  Liouville  operator  T(h)  accounts  for  the  spontaneous  decay.  Its 
definition  will  be  elaborated  in  the  next  section.  We  mention  that  T(h)  depends 
on  h,  which  embodies  the  second  geometry  dependence  of  the  equation  of  motion  for 
p(t). 

III.  SPONTANEOUS  DECAY 

Due  to  the  coupling  of  the  atomic  dipole  p  to  the  vacuum  of  the 
electromagnetic  field,  the  atom  decays  from  its  excited  state  to  its  ground 

29 

state,  together  with  the  emission  of  a  fluorescent  photon.  In  a  previous  paper 

we  evaluated  the  atomic  relaxation  operator  r(h)  for  an  arbitrary  atom  with 

degenerate  levels  and  near  a  perfect  conductor.  For  optical  transitions  we  can 

make  the  secular  approx imation^  (sometimes  called  'elimination  of  fast' 

31 

oscillating  terms'  )  and  neglect  the  imaginary  parts  of  the  relaxation 
constants.  If  we  subsequently  restrict  the  number  of  atomic  levels  to  two,  which 
are  coupled  by  the  laser,  then  T(h)  acquires  the  explicit  form 

r(h)o  -  i  £  ( 

m  m  m' 
e  g  e 

♦  Cm  .  m  <X)o|W<J.-;l> 

e  g  g  e 


•88® 


m'm'in  m 
e  g  g  e 


m  m  mm 
e  g  g  e 


L 

d  n _ 

e  g  e  g 


|j  m  ><J  m'|<j  m  lolj  m!>  , 

g  g  g  g*  e  «'  ®  « 


(3.1) 


which  defines  its  action  on  an  arbitrary  density  operator  o.  Through  the 
functions 


bx(x) 
bii  ( x) 


.  „fcos(2x)  sln(2x). 

(2x)2  (2x; 

.  3f sin(2x)  cos(2x)  sln(2x) , 

2‘  2*  <2x>2  ‘  (2x)2  '  ' 


(3.2) 

(3.3) 


the  relaxation  constants  can  be  expressed  in  terms  of  the  matrix  elements  of  the 
perpendicular  and  parallel  components  of  the  dipole  moment,  according  to 


U) 


c  ,  ,(h) 
m  m  mm 
e  g  g  e 


3irhe  c' 
o 


[  bx(  wh/  c  )  <  j  eme  I  yi  I  j  gng>  •  <  j  gmg  I  Hx  I J 

b|j  ( uoh/  c  )  <  J  e“e  I HJI I J  g°g>  *  <  J  gBi  I  H||  M  e®^ ]  * 


(3.4) 


The  presence  of  the  surface  is  fully  incorporated  in  the  functions  bx  and  b||, 
which  approach  unity  for  h  -*  •. 

Next  we  take  advantage  of  the  rotational  properties  of  the  atomic  wave 
functions  and  the  fact  that  p  is  a  vector  operator.  With  respect  to  the 
spherical  unit  vectors 


eA  -  e_  ,  •  -  ?(e_±ie  )//2  , 


:±1 


-x  -y- 


(3.5) 


we  can  decompose  y  ■  pj_  +  y..  as 


!!X‘ “on  •  H|-I  v!  • 


(3.6) 


T-±l 


7 


and  from 


t 

W 


p  it  follows  that 


t 


P 


T 


(-1)T 


w 


-T 


(3.7) 


Then  we  write  b^  ■  bj_  and  b±^  ■  by,  which  enables  us  to  cast  Eq.  (3.4)  in  the 

form 


u 


c  _  _,_.(h) 

m  ED  ID  m 

e  g  g  e 


3nhe  c' 
o 


\  bT(u,oh/c)<jemelWTlJg°g><Jg“gl\ljeme>  *  (3-8) 


Matrix  elements  of  p^  can  be  simplified  with  the  Wigner-Eckart  theorem,  which 
32 

states  that 


<JemeKIJgV  *  (JgmglTlJe“e)<je^Ml  I  Jg>//2Je+l  • 


8  8 


(3.9) 


and  the  matrix  elements  of  p^  then  follow  from  Eq.  (3.7).  Combining  Eqs.  (3.8) 
and  (3.9)  then  yields 


cm  m  »V(h)  "  AJ  bT(‘*>oh/c)(jgBglT,Je"e)(Jg,DilTljeBi)  * 

•  g  g  •  “  *  * 


(3.10) 


where  we  introduced  the  Einstein  coefficient  A^  for  a  free  atom  (no  surface) 


“o  l<J.IMIVl2 

*f  ‘  3,»c  c3  2V> 
o 


(3.11) 


Since  the  Clebsch-Gordan  coefficients  are  real,  all  relaxation  constants 

c  ,  ,  (h)  are  real.  Expression  (3.10)  reveals  the  surface  effect, 

■  ■mm 

e  g  g  e 

incorporated  in  the  parameter  functions  bT(o»oh/c)  on  the  decay  constants. 
Insertion  of  Eq.  (3.10)  into  Eq.  (3.1)  gives  T(h)  in  terms  of  Af,h  and 


geometrical  Clebsch-Gordan  coefficients.  The  coefficients  c  ,  ,(h)  are  only 

ID  D  B  n 

e  g  g  e 

non-zero  if  m  -a  •  m'-m',  as  follows  from  Eq.  (3.10),  so  for  m  ■  m'  we  have  m 
gege*  gg  e 

*  .  Therefore  the  summation  over  m'  in  the  first  two  terms  of  Eq.  (3.1)  is 

6  6 

trivial,  and  furthermore  we  notice  that  the  summation  over  involves  only  the 
decay  constants.  With  these  observations  the  expression  for  T(h)  reduces  to 

r(h)o  -  Am  (h)(|jeae><jeme|o  +  o|  j^Xj^J  )  -  Af£  b^uh/c) 


)  (j  o  It  I  j  m  )|j  m  ><j  m  lolj  m'XJ  m‘|(j  m'ltlj  m' )  ,  (3.12) 

L  g  8  e  e  ,Jg  g  Je  e1  ,Je  e  Jg  g1  Jg  g  ,Je  e 


m  m 
e  g 

m'm' 
e  g 


where  A  (h)  is  defined  by 


m 


A  (h)  -  aA  b  (u>  h/c)  J  (j  o  1t|J  m  )2  . 

m  f^to  L  g  g  e  e 


8 


(3.13) 


From  the  equation  of  motion  (2.8)  we  find  that  the  population  of  level 
b.V  obeys 


&  <Je»Jp(t)|jeae>  -  -<jeme|(r(h)p(t))|jeBe> 


(3.14) 


in  the  absence  of  a  laser  field,  and  from  Eq.  (3.12)  we  obtain 


-  Aa  (h)<J#m# j p(t) | j#*#> 


(3.15) 


Hence  the  population  of  Jj#m#>  decays  exponentially,  with  an  inverse  lifetime 

equal  to  A  (h),  which  elucidates  the  significance  of  this  parameter.  Without  a 

n 

e 

surface  we  have  J  1,  and  the  summations  in  Eq.  (3.13)  can  be  performed  to  give 
the  result 


A  (h)  ■  Af 
"e  f 


(3.16) 


Therefore,  the  population  of  an  atomic  level  |j  m  >  decays  in  free  space  with  A-, 

which  is  independent  of  m  .  This  in  turn  expresses  the  rotational  invariance  of 

spontaneous  decay  for  a  free  atom.  Conversely,  a  preferred  direction  in  space 

(the  z-axis),  which  is  imposed  by  the  presence  of  the  surface,  destroys  the 

rotational  symmetry,  which  implies  that  the  relaxation  depends  on  the  orientation 

of  the  atom  (the  m  -value  of  the  substate).  Prom  b  *  b  we  readily  derive 
e  t  -t 

K  <h>  “  A  m  00  .  (3.17) 

ID  “ID 

e  e 

which  expresses  the  remaining  rotational  symmetry  about  the  z-axis. 

IV.  TRANSFORMATION  AND  AVERAGE 

In  the  equation  of  motion  (2.10)  the  Hamiltonian  is  explicitly  time 
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dependent  through  the  external  field  E^(t).  With  the  standard  transformation 
of  the  density  operator 

-i(w.  t+$(t) )L 

p'(t)  •  e  8p(t)  ,  (4.1) 


where  L  is 
g 


the  commutator  with  the  projector  on  the  ground  state, 


e.g.. 


Lg°  *  *Pg’°*  *  (4.2) 

the  time  dependence  is  greatly  simplified.  We  find  the  equation  for  p'(t)  to  be 

•  <Ld<h)  ♦  Lg  -  ir(h»p*  (4.3) 

In  the  compact  Liouville  notation.  Oscillations  with  the  optical  frequency 
have  disappeared,  and  the  phase  fluctuations  only  enter  as  a  multiplicative 
factor  d^/dt  to  p'(t),  which  is  the  standard  form  of  a  stochastic  differential 
equation.  Here  the  Liouvillian  L.(h)  is  defined  as 


Ld(h)o  -  #  1(Hd(h),o] 


(A. 4) 


in  terms  of  the  dressed-atom  Hamiltonian 


Hd(h)  -  Ha  +  ftuijP  -  |h(n(h)  +  n(h)  +  )  , 


(4.5) 


which  has  the  significance  of  the  atom,  dressed  with  the  photons  of  the  laser  and 

including  the  dipole  interaction.^ 

Due  to  the  term  L  d$/dt  in  Eq.  (4.3),  the  density  operator  p'(t)  is  a 
8 

stochastic  process.  Only  its  average 


n(t)  -  <p'(t)> 


(4.6) 


has  relevance,  and  the  equation  for  Il(t)  reads' 


ldt  "  (Ld(h)  '  iW  "  ir(h))n  » 


(4.7) 


with 


W  -  AL 


(4.8) 


in  terms  of  the  half-width  at  half  maximum  A  of  the  Lorentzian  laser  profile. 

Both  the  laser  linewidth  (W)  and  spontaneous  decay  (r(h))  give  rise  to 
dealing  of  the  free  evolution  (Ld(h))  of  the  dressed  atom.  After  a  transient 
time  of  the  order  of  (A  +  A^)  *,  following  the  switch-on  of  the  laser,  the  ate 
will  have  reached  a  steady  state.  If  we  write 


I  -  lia  n(t)  , 

t-»- 


(4.9) 


then  this  H  is  the  solution  of 


(Lj(h)  -  iW  -  ir(h))II  -  0  . 


(4.10) 


a*  r  *  ) 


Since  the  external  field  is  assumed  to  be  a  CW  laser,  we  will  always  consider  the 
atom  to  be  in  the  stationary  state  1!. 


V.  FLUORESCENCE 

Continuous  excitation  of  the  upper  states  |j  m  >  by  the  laser  and  subsequent 

spontaneous  decay  to  the  ground  states  |j  m  >  goes  together  with  the  emission  of 

8  8 

a  stationary  (fluorescence)  radiation  field.  This  signal  is  detected  by  a 
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photomultiplier  tube  (PM)  at  a  distance  r  from  the  atom.  In  a  previous  paper 
we  derived  the  general  expression  for  the  fluorescence  in  the  far  field,  which 
reads 


Ef(r,t)  «  - r~ 

4ire  c  r 
o 


(r*(p(t+T-r/c)  +  u' (t-T-r/c))]  *  r  . 


(5.1) 


Here,  r  *  r/r  is  the  direction  of  propagation,  and  t  ■  r-e  h/c  equals  half  the 
delay  time  between  the  directly-emitted  wave  and  the  wave  which  is  first 
reflected  by  the  surface.  The  time  dependence  of  the  operator  p  represents  the 
Heisenberg  picture  according  to 


p(t)  »  e 


iL(t-t  ) 


(5.2) 


for  any  time  interval  t  -  tQ  and  in  terms  of  Lo  ■  Ji  1(H,o],  with  H  tha 
Hamiltonian  of  the  entire  system.  If  we  divide  the  operator  p  as  p  »  p^  ♦  yg, 
then  p1  is  defined  as 


y’  "  el  -  y« 


(5.3) 


VI.  DETECTION 


From  Eq.  (5.1)  it  follows  that  the  polarization  of  the  fluorescence  is 


12 


determined  by  the  vector  operators  y  and  y',  which  have  a  different  direction  in 

space  and  a  different  retardation  time.  In  order  to  keep  track  of  the 

polarization  properties  of  the  radiation,  we  suppose  that  the  radiation  passes  a 

polarizer,  before  entering  the  PM,  which  transmits  radiation  with  polarization 
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€p.  Prom  the  theory  of  photon  detection  by  a  PM,  '  we  then  know  that  the 
detector  responds  to  the  positive  frequency  part  of  the  scalar  quantity 


£(t)  -  Ef ( r , t ) * e*  , 


(6.1) 


which  is  the  component  of  the  fluorescence.  Then  we  position  the  polarizer  in 

such  a  way  that  the  angle  of  incidence  is  zero  (e.g.,  r-e^  *  0),  and  we  suppress 
the  overall  retardation  with  r/c.  We  then  find 


I2  1Lt  4f 

£(t)  -  — ~ — y  [e1  T  y  + 

Aire  c  r 
o 


-ILt  . 

e  h  1**5 


(6.2) 


where  y  and  y'  are  time  independent  operators  (Schrodinger  picture).  Then  we 

notice  that  the  surface  can  only  have  a  significant  effect  if  the  distance 

between  the  atom  and  the  xy-plane  is  a  few  optical  wavelengths.  Hence  the 

retardation  t  in  Bq.  (6.2)  is  of  the  order  of  a  few  optical  cycles  (t  -  u  *1), 

o 

which  implies  that  ve  can  neglect  the  interaction  term  with  the  laser  in  L  in  the 

exponentials.  In  physical  terms  this  means  that  the  laser  cannot  cause  a 

transition  between  the  upper  and  lower  states  on  a  time  scale  t,  which  is 

2 

completely  justified  for  optical  lasers.  Furthermore,  the  L  is  merely  an 

overall  factor,  which  hardly  affects  the  time  evolution  of  £(t).  Therefore,  we 

2 

can  safely  make  the  substitution  Lo  ■*  [H^.oJ/h  in  expCtiLr)  and  L  .  Combining 
everything  then  yields  for  the  positive  frequency  part  of  S(t) 


(6.3) 


<£(*■) 


(jJ 

(t)  *  - ^“5“  eiLt{P  (u  (o>  t)*-C*)P  } 

4 ire  c2r  8  ~s  0  *P  6 

o 


VII.  PHOTON  CORRELATIONS 

If  the  active  part  of  the  PM  is  exposed  to  electromagnetic  radiation  it  will 
absorb  photons  from  this  incident  field.  Time  correlations  between  various 
photons  are  most  conveniently  expressed  in  the  function  I^t^.t^dtjdt^  which 
has,  by  definition,  the  significance  of  the  probability  for  the  detection  of  a 
photon  in  [tj,t^+dt^)  together  with  a  detection  in  [t2*t2+dt2],  and  irrespective 
of  detections  at  other  times.  The  basic  relation  between  and  the  incident 
radiation  is  given  by^’^ 


c2<^(*)(t1)^(')(t2)^(+)(t2)^(+)(t1» 


(7.1) 


for  t2  >  tj.  Here,  C  is  a  detector  parameter  (depending  on  efficiency,  aperture, 
etc.)  and  the  angle  brackets  indicate  a  quantum  average.  For  our  case,  where  the 
laser  field  is  taken  as  a  classical  wave  with  a  stochastic  phase,  the  brackets 
also  include  an  average  over  the  process  ♦(t).  The  quantity  <...>  in  Eq.  (7.1) 
is  usually  referred  to  as  the  intensity  correlation  of  the  field.  For  quantum 
radiation,  like  fluorescence,  the  different  factors  in  brackets  do  not  commute, 
which  gives  rise  to  non-trivial  correlations  between  the  detections  of  photons. 

It  is  a  standard  procedure  to  work  out  the  intensity  correlation  and  to 
express  the  time  regression  from  t^  to  t^  in  the  time-evolution  operator  for  the 
density  operator 


j  U(t) 

t 


■  e 


» 


-i(Ld(h)-iV-ir(h))t 


(7.2) 


and  the  photon-emission  operator  R,  defined  by 
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Ro  -  5*+*(0)o^“*(0)  .  (7.3) 

If  we  assume  that  the  atom  has  reached  its  steady  state  H,  then  the  intensity 
correlation  function  is  found  to  be 

I2(tl'tl+T)  “  C2TrRU(-r)R5  -  f(t)I  ,  x  >  0  ,  (7.4) 

which  defines  the  function  f(x).  Here,  I  equals  the  uncorrelated  intensity  of 
the  radiation,  which  can  be  expressed  as 

I  -  CTrRH  .  (7.5) 

For  a  zero  delay  time  x  we  have  U(0)  ■  1  and  f(0)  »  C2TrR2H/I.  With  Eq. 
(6.3)  we  obtain 

R2  -  0  ,  (7.6) 

and  therefor* 

f(0)  -  0  .  (7.7) 

This  famous  relation  for  the  detection  of  two  fluorescent  photons  is  termed 
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antibunching.  In  the  case  that  the  delay  time  x  is  very  long,  we  find 


f(-)  -  I  . 


(7.8) 


which  shows  that  the  second  detection  is  independent  of  the  first  one,  e.g.,  the 
detections  are  uncorrelated.  Equation  (7.8)  follows  from  the  identity 


lim  U(t)o  ■  IlTro  (7.9) 

t-*» 


for  any  o,  which  states  that  every  o  evolves  to  a  unique  steady  state  II  and  that 
the  trace  of  o  is  conserved  in  the  time  evolution  with  U(t). 

With  Eqs.  (7.3),  (6.3),  (2.3)  and  (2.4)  we  find  for  the  emission  operator 
the  general  expression 

Ro  *  )  |j  m  ><J  m1 I  <j  m  lolj  m'> 

L  8  8  8  8  e  e'  e  * 

m  m 
e  g 

minig 


< j  m  | u  (u>  t  )  •  e  |  j  m  >*  < j  m'  J u  ( w  t  )  •  e  |  j  ■*  > 
e  e'-s  o  -p'  g  g  e  e'-s  o  -p 1  g  g 


(7.10) 


where  an  overall  factor 
detector  parameter  (. 
with  Eq.  (3.9). 


2  2  2 

vu  /4ne  c  r)  is  omitted,  since  it  can  be  absored  in  the 
o  o 

Matrix  elements  of  p  (t^i)  can  be  evaluated  explicitly 


VIII.  1  ■  0  ->  1^  ■  1  TRANSITION 

In  order  to  study  the  effect  of  the  presence  of  the  metal  surface  on  the 

photon  correlations,  we  have  to  prescribe  more  details  of  the  system.  In  the 

remainder  of  this  paper  we  consider  the  important  example  J  »  0  ■  1.  For 

j  -  0  we  have  a  ■  0  only,  and  in  view  of  Eq.  (3.9)  this  non-degenerate  state 
8  8 

couples  only  with  J#  ■  1,  which  is  three-fold  degenerate.  We  introduce  the 


abbreviations 


(8.1) 


|jg®g>  "  |0  0>  -►  |g>  , 

|j  m  >  ™  1 1  m  >  -►  | *  n  “  ®  ■  -1.0.1  .  (8.2) 

6  6  6  6 

Spontaneous  decay  of  this  atom  is  described  by  the  Liouvillian  T(h)  from  Eq. 
(3.12),  and  the  inverse  lifetimes  of  the  states  |n>  are  given  by  Eq.  (3.13).  The 
occurring  Clebsch-Gordan  coefficients  all  have  the  form 

(0  0  1  x 1 1  n)  -  6  ,  (8.3) 

which  yields  for  the  relaxation  constants 

A  (h)  -  A,b  (u>  h/c)  .  (8.4) 

n  f  n  o 

We  observe  that  the  substates  |n  ■  *1>  decay  with  AjbijfuMi/c),  which  equals  the 
inverse  lifetime  of  a  dipole,  oriented  parallel  to  the  xy-plane.  The  substate 
|n  *  0>  has  the  Einstein  coefficient  A^b^c^h/c) ,  and  therefore  it  decays  as  a 
perpendicular  dipole.  Notice  that  the  differences  in  decay  constants  of  the 
various  substates  are  merely  determined  by  the  geometry  of  the  system.  A  two- 
state  model-atom  treatment  of  fluorescence  near  a  surface  leaves  the  dipole 
direction,  and  thereby  the  lifetime  of  the  excited  state,  undetermined.  One  can 
then  choose  the  direction  of  the  matrix  element  <e|p|g>  arbitrarily,  and 
subsequently  average  results  over  a  randomization  of  this  direction.  It  should 
be  obvious  now  that  this  is  not  a  correct  procedure. 

The  spontaneous -decay  operator  for  this  system  attains  the  simple  form 

r(h)o  ■  J^AJi(h)(|n><n|o  ♦  ojn><n|  -  2|g><g|<n|o|n>)  »  (8.5) 

n 

and  it  involves  only  the  Einstein  coefficients  A  (h)  for  the  three  substates. 

n 


IX.  ATOMIC  STATE 


In  this  section  we  expand  the  equation  of  motion  (4.7)  for  the  atomic  state 
0(t)  in  matrix  elements.  The  external  field  enters  through  the  Rabi  operator 
Q(h),  Eq.  (2.9),  and  we  shall  abbreviate  its  non-zero  matrix  elements  as 


n^h)  -  <n  | n(h)  | g>  , 

(9.1) 

which  are  explicitly 

5  E 

no(h)  ■  75  ~i  <1ll^ll°>  co*(hk,)!-;0  * 

(9.2) 

-  E 

n±(h)  ■  73  <1 1 1  y 1 1 0>i  sin(hk2)e-e*1  . 

(9.3) 

With  the  definitions  of  L.(h),  W  and  T(h)  we  then  obtain 

<n | n | n ' >  -  -  i(An  +  An« )<n|n|n*>  *  ii(n*,<n|n|g>  -  nn<g|n|r1'>)  ,  (9.4) 

<g|n|g>  -  £  An<n|n|n>  -  ii£(on<g|n|n>  -  n*<n|n|g>)  ,  (9.5) 

n  n 

<n | n | g>  -  -(iAn  +  a  -  iA)<njn|g>  +  ii(nn<g|n|g>  -  £  nn,<n|n|n*>)  ,  (9.6) 

n' 

<g|n|n>  •  -(jAn  ♦  x  ♦  it)  <g|n|n>  -  Ji(Q*<g|n|g>  -  £  Q*,<n'|n|n»  .  (9.7) 

n' 

where  the  dependence  of  and  0^  on  h  is  suppressed,  and  we  Introduced  the 
detuning  of  the  laser  from  resonance  by 


(9.8) 


If  we  make  all  the  combinations  with  q  and  n ' ,  then  Eqs.  (9.4)-(9.7)  constitute  a 
set  of  16  coupled  first-order  linear  differential  equations. 
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X.  INTENSITY  CORRELATIONS 

For  the  J  ■  0  -*■  j  ■  1  transition,  the  photon-emission  operator  R,  Eq. 

s  ® 

(7.10),  reduces  to 

R°  "  P8  I  <r» 1 0 1 1  * ><n  |  jus(u)0T )  -  |8>*<n '  I  ys(a)oT )  * ep  I g>  ,  (10.1) 

nn' 

and  the  matrix  elements  of  y  (w  x)  are 

~s  o 


<n|Ms(u,oT)|g> 


[cos(u,oT)e05no  +  i 


sin(u)  t)" 
o 


l 


;*.«nn.)  <1| |p| |0>  2//3  . 

(10.2) 


We  notice  that  the  summation  over  n  and  n'  in  Eq.  (10.1)  is  Just  a  number  (not  an 
operator),  and  therefore  Ro  is  proportional  to  the  projector  on  the  ground 
state  for  every  o.  This  reflects  the  fact  that  after  every  emission  the  atom  is 
necessarily  in  its  ground  state  |g>.  Since  for  Jg  “  0  this  state  is  non¬ 
degenerate,  the  final  state  after  a  photon  emission  is  unique,  or  in  other  words, 
independent  of  the  history  of  the  atom.  If  the  ground  state  were  degenerate, 
then  we  observe  from  Eq.  (7.10)  that  the  density  operator  of  the  atom  after 

emission  is  a  superposition  of  (Liouville)  substates  |j  m  ><J  m' I ,  where  the 

8  8  8  8 

distribution  over  states  depends  on  the  matrix  elements  of  o,  which  equals  the 
density  operator  Just  before  the  emission. 

Combination  of  Eqs.  (7.S)  and  (10.1)  gives  for  the  intensity  of  the 
fluorescence 

1  "  c  I  nnl<n^s(V)*5p*8>l2  »  (10.3) 

n 

in  terms  of  the  steady-state  populations  n^  ■  <n|H|n>  of  the  upper  levels.  The 
tvo-photon  correlation  f(t)  then  becomes 

8  ’ 


f(t)  -  CTrRU(t)P 


(10.4) 


which  involves  only  a  single  photon-emission  operator  R,  rather  than  two  as  in 
Eq.  (7.4).  This  is  again  due  to  the  fact  that  the  ground  state  is  non¬ 
degenerate.  After  the  first  emission  the  density  operator  equals  P  ,  which  is 

8 

the  right-most  factor  in  Eq.  10.4). 

Substitution  of  the  exponential  for  U(t),  Eq.  (7.2),  into  Eq.  (10.4)  gives 
an  explicit  expression  for  the  correlation  function  f(t).  The  Liouvillian 

L.(h)  -  iW  -  ir(h),  however,  is  a  16  *  16  matrix,  which  makes  the  evaluation  of 

a 

the  exponential  intractable,  in  general.  In  most  practical  cases  (see  the  next 
section),  the  computation  of  f(t)  is  facilitated  if  we  first  transform  to  the 
Laplace  domain.  With 

?(s)  ■  f  dt  e  Stf(t)  ,  (10.6) 

J0 

we  find  immediately 

f(s)  .  CTrR  - lC^  r>r;  m  Pg  .  (10.7) 

which  only  requires  the  inversion  of  the  matrix.  Of  course,  the  result  must  be 
transformed  back  to  the  time  domain  in  order  to  achieve  an  explicit  expression 
for  f(t). 

XI.  POLARIZATION 

To  fix  the  geometry  of  the  experiment,  there  are  four  directions  which 

remain  to  be  specified.  These  ares  (1)  the  laser  polarization  c,  (2)  the  laser 

propagation  direction  (  (restricted  by  £*c  ■  0),  (3)  the  observation  direction 

and  (4)  the  polarization  unit-vector  e  (restricted  by  r*c  ■  0).  We  shall 

-P  ~  ~P 

consider  two  complementary  situations,  for  which  we  evaluate  the  two-photon 
correlation  function  explicitly. 


(a)  Linear  polarization.  The  laser  polarization  is  chosen  as  e  *  e^  ■  e  ,  for 
which  the  direction  of  propagation  fc  must  be  in  the  xy-plane.  Then  we  have 

k  ■  0,  and  from  Eqs.  (9.2)  and  (9.3)  we  find  that  the  only  non-zero  matrix 
element  of  the  Rabi  operator  is 

-  E 

°o(h)  "  75  <JI MJ0>  *  (n.i) 

which  is  independent  of  the  distance  h  between  the  atom  and  the  surface.  Prom 
Eqs.  (9.4)-(9.7)  it  follows  that  this  external  field  couples  only  the  ground 
state  |g>  with  the  upper  state  |0>, which  has  the  decay  constant  Ag(h)  ■ 
Ajb^(u>oh/c) .  The  PM  will  be  positioned  in  such  a  place  that  it  detects 
radiation,  which  is  emitted  along  the  surface,  e.g.,  r  is  in  the  xy-plane.  Then 
we  filter  the  fluorescence  with  *  e^.  For  this  configuration  the  photon- 
emission  operator  becomes 

RQo  -  Pg  <0 | a | 0>  ,  (11.2) 

where  a  factor  |<1 j |u| |Q>|24/3  is  absorbed  in  C. 

(b)  Circular  polarization.  In  this  case  we  take  e  ■  e  ,  ■  e  ,  i  »  -eft  ■  -r  and 

-  ~il  -p  -0 

therefore  k  »  -k  -  -w  /c.  The  Interaction  with  the  laser  field  now  becomes 
z  —  o 

2  * 

0±1(h)  ■  75  <1||m||0>  i  sin(-hk)  ,  (11.3) 

which  couples  only  |g>  with  |il>.  Hence  the  relevant  Einstein  coefficient  equals 
A±^(h)  ■  A^b^(uoh/c),  and  for  the  operator  R  we  obtain 

R±1o  -  Pg  <±1 1 o |  ±1>  sin2(u»oh/c)  .  (11.4) 

If  we  denote  the  wavelength  of  the  radiation  by  A  ■  2w/k  -  2vc/w  ,  then  we  notice 

—  o 

that  both  0^(h)  and  R^  vanish  for  h  ■  n  *  A/2,  with  n  ■  0,1,2,...  .  For 
circular  polarization  the  fluorescence  intensity  is  identically  zero  if  h  equals 


*rV 


an  integer  multiple  of  A.  Furthermore,  the  emission  operator,  the  Einstein 
coefficient  and  the  Rabi  frequency  are  identical  for  left-  and  right-handed 
polarizations,  and  so  the  photon  correlations  are  identical  (provided  of  course 


2 


that  we  take  e  ■  c^),  which  expresses  again  the  rotational  symmetry  about  the  z- 
axis. 

Now  it  is  an  easy  matter  to  evaluate  f(s),  Eq.  (10.7),  for  both  linear  and 

circular  polarization  simultaneously.  The  set  of  equations  (9.4)-(9.7)  defines 

the  matrix  s  +  iL.(h)  +  W  +  r(h),  and  with  some  algebra  we  then  find 
a 


ms+m  +A)|nnr 


f  (s)  =  i  *  - n - u - 

n  S  in  I 2( s+jA  +A)  +  (A  +s)[(s+jA  +A)2  +  A2] 

n  q  q  q 


[ 1  -  6  .  cos  (u  h/c))  , 

q.tl  o 


(11.5) 


for  q  ■  -1,0,1.  Then  the  fluorescence  intensity  follows  from 


I  ■  lim  f  (t)  •  lim  sf  (s)  . 


(11.6) 


XII.  SPECIAL  CASES 

Transformation  of  the  general  result  (11.5)  for  ?  (s)  to  the  time  domain  is 

h 

straightforward,  but  the  expression  for  f^(t)  is  very  untransparent.  A  case  of 
practical  Importance  is  the  low-intensity  limit,  where  the  coupling  strength 


between  the  dipole  and  the  laser  field 


^  *  7§  \  !<MMI°>I 


(12.1) 


tends  to  zero  in  comparison  with  the  relaxation  parameter  A^.  Normalized  with 

the  uncorrelated  intensity  I  ■  f  («•),  f  (t)  is  found  to  be 

Oh  h 


X  -  |A  (A+U  )2  +  A2 


f  (t)/I  -  1  -  -  -  r-'1  - °-r - J 

q  q  X  +  JA  (X.|A  )2  +  a2 

h 


e 


1 


e 


+ 


A 

_ n 

X  +  JA 

n 


(A-iAn)2  +  A2 


-(A+jA^t 


*  ((A2  -  $A2  -  A2 )cos( At)  -  2AAsin(At)] 

n 


(12.2) 


Figures  1  to  4  illustrate  the  dependence  of  f  (t)/I  on  the  parameters  A  and  A, 
the  laser  polarization  and  the  distance  of  the  atom  to  the  surface. 

For  higher  irradiances  we  assume  that  is  in  very  close  resonance  with  u»o, 
e.g.,  we  take  A  »  0.  Then  we  obtain 


f  (t)/I  »  1  -  e  n  ( cos( 8  t)  +  r11  sin(B  t)) 

n  n  n  Pn  n 


(12.3) 


with 


-  3A^/ 4  +  iA  ,  (12.4) 

pn  *  -  i(Un  -  A)2]1  ,  (12.5) 

and  typical  behavior  is  illustrated  in  Figs.  5  and  6. 

XIII.  CONCLUSIONS 

An  excited  atom  near  a  surface  has  a  different  lifetime  1/A  than  in  free 
space.  Standard  methods  to  measure  A  and  its  dependence  on  the  distance  of  the 
•tom  to  the  surface  are:  (1)  A  equals  the  full  width  of  the  low  intensity 
absorption  profile  and  (2)  l/A  equals  the  delay  time  in  the  emission  of  a 
fluorescent  photon  after  excitation  of  the  atom.  In  this  paper  we  considered  the 
situation  where  the  atom  is  continuously  driven  by  a  laser  field,  and  we 
evaluated  the  temporal  correlation  between  photons  of  the  emitted  fluorescence. 

Alter  presentation  of  the  general  theory,  we  focussed  on  the  j  -  0  J  *1 

^  8  • 

transition.  In  the  most  simple  case,  where  the  irradiance  is  monochromatic  (A  ■ 
0),  exactly  on  resonance  (A  «  0),  sufficiently  weak  (n^  0), 


and  has  a  polarization  n>  we  found  (Eq.  (12.2)) 


-*A  t  o 

f  (t)  -  I  (1  -  e  n  )2  .  (13.1) 

n  n 

We  conclude  that  the  two-photon  correlation  function  f(t),  which  is  accessible  to 
observation  by  a  photon  counter,  is  essentially  determined  by  the  lifetime  1/A^ 
of  the  driven  upper  substate  of  the  atom.  This  result  will  possibly  provide  an 
alternative  optical  technique  for  the  investigation  of  atoms  near  a  metal 
surface. 
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FIGURE  CAPTIONS 

j?  L&-! _ I  •  Plot  of  the  two-photon  correlation  function  f^(t)/I^  in  the  low- 

intensity  limit.  Time  and  frequencies  are  given  in  units  of  1/A^  and  A^, 
respectively,  and  the  distance  between  the  atom  and  the  surface  h  is  in  units  of 
an  optical  wavelength  2irc/wo.  In  this  figure  we  take  h  ■  0.2,  A  ■  0.5  and  linear 
polarization.  Curves  a,  b  and  c  illustrate  the  dependence  on  the  detuning,  with 
A  *  0,  2  and  4,  respectively.  Notice  that  f(t)  can  become  larger  than  f(»), 
which  implies  an  enhancement  for  the  detection  of  a  photon  at  time  t  after  the 
detection  of  a  photon  at  time  zero,  in  comparison  with  the  uncorrelated  detection 
probability  I. 

Fig.  2.  Same  as  Fig.  1  but  now  we  vary  the  laser  linewidth  as  A  *  0.2,  0.6  and 
3.0  for  curves  a,  b  and  c,  respectively,  and  the  detuning  is  fixed  at  A  ■  4.  For 
a  large  bandwidth,  curve  c,  the  oscillatory  structure  is  washed  out  entirely. 

Low-intensity  f^(t)  for  h  ■  0.2,  A  ■  0.1  and  A  ■  2.  Curves  a  and  b 
correspond  to  linear  and  circular  laser  polarization,  respectively. 

Fig.  4.  Two-photon  correlation  in  the  low-intensity  limit  with  A  ■  0.2,  A  »  2 
and  circular  polarization.  This  picture  illustrates  the  dependence  on  the  atom- 
surface  distance  h,  with  h  ■  0.05  and  0.1  for  curves  a  and  b,  respectively.  In 
curve  a  the  approach  to  f^(«)/I^  ■  1  for  Increasing  time  occurs  very  slowly. 

This  can  be  understood  from  the  fact  that  for  circular  polarization  the  atom 
behaves  effectively  as  a  parallel  dipole.  Its  matrix  element,  which  determines 
the  inverse  lifetime,  tends  to  zero  for  h  ■+  0,  and  therefore  the  lifetime  itself 
becomes  very  large.  This  in  turn  prohibits  the  emission  of  the  second  photon  at 
t  >  0,  after  emission  of  the  first  photon  and  subsequent  excitation  by  the  laser. 

Fig.  5.  Correlation  function  for  zero  detuning,  h  ■  0.1,  A  ■  0.5  and  ■  3. 
Curves  a  and  b  correspond  to  linear  and  circular  polarization,  respectively. 


2 


Fig.  6.  The  parameters  in  this  figure  are  X  ■  0.5,  A  ■  0,  0^  ■  8,  and  the  laser 
is  circularly  polarized.  Curves  a  and  b  represent  the  change  in  correlation  as  a 
function  of  h.  Vie  take  h  *  0.1  and  1  for  a  and  b  respectively.  It  appears  that 
the  atom-surface  distance  can  affect  the  photon-correlation  function 
dramatically,  which  implies  that  a  measurement  of  f^(t)/I^  might  appear  to  be  a 
sensitive  technique  for  the  determination  of  h  of  atomic  adsorbates. 
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